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Abstract 

We present a combined analysis of XMM-Newton, Chan- 
dra and Rosat observations of the isolated neutron star 
RX J0720.4-3f 25, spanning a total period of ~ 7 years. 
We develop a maximum likelihood periodogramme based 
on AC-statistic and maximum likelihood method, which 
are appropriate for sparse event lists. As an a posteriori 
check, we have folded a further BeppoSAX-dataset with 
the period predicted at the time of that observation, find- 
ing that the phase is consistent. 

The study of the spin history and the measure of the 
spin-down rate are of extreme importance in understand- 
ing the mechanism powering RX J0720. 4-3125. The value 
of P, here measured for the first time, is « 10 -14 s/s and 
can not be explained in terms of propeller or torque from 
fossil, disk, When interpreted in terms of dipolar losses, 



it gives a magnetic field of B s-i 10 13 G, making also im- 



Until a few years ago, dim NSs were thought to con- 
stitute a class stand alone and two mechanisms were pro- 
posed for their emission: either accretion from the inter- 
stellar medium onto an old NS or release of thermal radi- 
ation from a younger, cooling object. More recently, based 
on the similarity of the periods, it has been suggested 
a possible evolutionary link between dim NSs, anoma- 
lous X-ray pulsars (AXPs), and soft gamma-ray repeaters 
(SGRs). Two kind of "unified" scenarios have been then 
proposed. In the first one, the three classes of objects are 
powered by dissipation of a decaying, superstrong mag- 
netic field (B > 10 14 - 10 15 G). In this case dim NSs are 
the descendants of SGRs and AXPs, and RX J0720.4-3125 
may be one of the closest old-magnetars. Alternatively, all 
the three classes may contain NSs with lower (canoni cal) 
magnetic field ( B « 10 12 G) endowed by a fossil disk (Al- 
par et al. 2001). In this case dim NSs in the propellor 



plausible that the source is accreting from the underdense 
surroundings. We also find unlikely that the field decayed 
from a much larger value (B w 10 15 G in the past, as 
expected in the "old magnetar" interpretation) since this 
scenario predicts a source age of ~ 10 4 yrs, too young to 
match the (low) observed X-ray luminosity. The observed 
properties are more compatible with a scenario in which 
the source is w 10 6 yrs old, and its magnetic field has not 
changed substantially over the lifetime. 

Key words: Stars: neutron — stars: oscillations — pulsars: 
general — magnetic fields. 



1. Introduction 



phase would be the progenitors of AXPs and SGRs, the 
latter having entered an accretion phase. 



Recently, |Paerels et al. (2001)| presented XMM-Newton 
spectra of RX J0720. 4-3125. The absence of electron or 
proton cyclotron resonances in the RGS range excluded 
magnetic fie lds of B w (0.3- 2) x 10 11 G and (0.5 - 2) x 



10 14 G (see ftane et al. 200 ip. Based on the same XMM- 



Newton observation, |Croppcr et al. (200 1)| presented the 
pulse-shape analysis. They derived an upper limit on the 
polar cap size, showing that an emitting region larger than 
~ 60° — 65° can be rejected at a confidence level of 90%. 
Whatever the mechanism, the X-ray emitting region is 
therefore confined to a relatively small fraction of the star 
surface. They also found that the hardness ratio is soft- 
est around the flux maximum. The same has been later 



RX JQI 720.4 3125 m a nearby, isolated neutron sta r (NS) 
detected by ROSAT during a Galactic plane survey (Haberl 
et al. 1997) and recently re-observed with XMM-Newton 
on 2000 May 13 QPacrcls et al. 200l| , proppcr et al. 200l| ) 
and 2001 November 21. The source exhibits all the com- 
mon characteristics of the other six ROSAT NS candi- 



discovered by Perna et al. (2001) in some AXPs. Cropper 



et al. (2001) suggested two possible explanations for this 



dates (hereafter dim NSs, see Treves et al. 200C for a re- 
view): a blackbody-like spectrum with kT ~ 80 eV; a 
large X-ray to optical flux ratio; a low X-ray luminosity, 
Lx ~ 10 30 — 10 31 erg/s; a low column density and no ev- 
idence for a binary companion. In addition, RX J0720.4- 
3125 is pulsating with a period P ~ 8.4 s. 



effect: either radiation beaming (as in their best-fitting 
model) or the presence of a spatially variable absorbing 
matter, co-rotating in the magnetosphere. The latter may 
be indeed the case if the star is propelling matter outward 
( [Alpar et al. 200 1| ). 

Further information about the nature of this puzzling 
source can be obtained by the spin history. Magnetars will 
spin-down at a rate P « 10" n (B/10 14 G) 2 /B ss" 1 , due 
to magneto-di polar losses. The prel iminary measure of P 
published by |Haberl et al. (1997)| for RX J0720.4-3125 
is uncertain to a considerable large value, and does not 



Proc. Symposium 'New Visions of the X-ray Universe in the XMM-Newton and Chandra Era' 26-30 November 2001, 



2 



even allow spin-up and spin-down to be discriminated. An 
accurate determination of P is therefore crucial, as well a 
tracking of the spin history of the source. Here we present 
a combined analysis of XMM-Newton, Chandra and Rosat 
data, spanning a period of ~ 7 years. For all details we 



refer to the paper Zane et al. (2002) 



2. Timing analysis 

The different observations used in our analysis are shown 
in table [|; the major datasets are from the two XMM ob- 
servations and from the 1996 Nov. 3 Rosat pointing, while 
the 1998 Rosat and the Chandra observations are valuable 
by nature of their several day durations. Our data origi- 
nate from instrumentation with widely differing sensitivi- 
ties: typical count rates vary from 1 count every ~ 3 s for 
Rosat HRI to ~ 6 counts/s for XMM-Newton PN. How- 
ever, none of these count rates is sufficiently high for a 
normal distribution of counts to be expected, thus stan- 
dard discrete Fourier Transforms are not directly applica- 
ble. For sparse data and event list data, we used instead 
Rayleigh Transform (i.e. de Jager 1991, Mardia 1972). It 
is also crucial for us to define precisely the confidence in- 
tervals to the derived quantities, in particular the period 
P. We do this by constructing MLP (maximum likelihood 
periodogrammes) which make no assump tions on da ta dis- 
tribution, and using the AC-statistics ( ]Cash 1979 ). The 
uncertainty in the period and the \ 2 can be read directly 
from the y-axis of the MLP (see figure @). 



Table 1. The ROSAT, Chandra and XMM-Newton ob- 
servations of RX J0720. 4-3125 used in this analysis. The 
entry in the third column is the effective exposure. 



Date 


Instrument 


Eff. Exp. (s) 


Label 


1993 Sep 27 


Rosat PSPC 


3221 


R93 


1996 Apr 25 


Rosat HRI 


3566 


R96a 


1996 May 7 


Rosat HRI 


3125 


R96b 


1996 Sep 27 


Rosat HRI 


1409 


R96c 


1996 Nov 3 


Rosat HRI 


33569 


R96d 


1998 Sep 27 


Rosat HRI 


3566 


R98 


2000 Feb 1 


Chandra HRC-S 


37635 


ChOO 


2000 May 13 


XMM MOS1 


61648 


XOOa 




MOS2 


61648 






PN 


62425 




2000 Nov 21 


XMM MOS1 


17997 


XOOb 




MOS2 


17994 






PN 


25651 




1997 Mar 16 


SAX LECS 




S97 



We first performed an MLP assuming P = on each 
of the longer pointing: R93, R96d, XOOa, XOOb (figure 1). 
There is no ambiguity in the period determinations and 
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Figure 1. (Left) Maximum likelihood periodogrammes 
(MLP) for three long datasets, R96d, XOOa (PN) and 
XOOb (PN), showing the periodicity at 8.391 s. These con- 
strain the selection of the strongest and second- strong est 
dips in the MLPs for the R98 and ChOO datasets re- 
spectively (right). The vertical line denotes a period of 
8.39113 s. The 68% and 90% confidence levels are at 



X' 



1.0 and 2.71 for one degree of freedom. 




Figure 2. The 68%, 90% and 99% contours for a linear 
least squares fit of R93, R96d, R98 XOOa, XOOb and CHOO 
(continuous elliptical regions). Parallel lines are the 68% 
and 90% contours of XOOa PN; tiny elliptical regions are 
the 68% and 90% contours for the combined R93 and R96 
datasets (see zoom). 



a linear least square fit using the 68% formal errors in 
the MLP gives P = 8.39113 ± 0.00011 s, P = 0.0 ± 5.5 x 
10 -13 s/s (here and in the following Pq is referenced to the 
start of the R93 run). This upper limit on P permits an 
unambiguous determination of the peaks in the ChOO and 
R98 power spectra. Adding these to the linear square fit 
gives P = 8.39107 ± 0.00005 and P = 2.7 x 10~ 13 ± 2.5 x 
10" 13 . The 68, 90 and 99% confidence intervals are shown 
in figure @, as well the 68 and 90% intervals derived from 
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XOOa. With the improved (Pq,P) values, we performed 
an MLP on the combined R93 and R96 datasets. As a 
result, the confidence contour break up into small region 
(aliases) in the (Pq , P) plane (see zoom in figure 2). With 
this further restriction, we finally do the MLP on all data. 

We derive two pairs of values (Pq , P) which cannot 
be further discriminated between on statistical grounds 
(table g). 

Table 2. The two pairs of best fitting (Pq , P) values. A\ 2 
is the difference between the \ 2 of solution (2) and that of 
solution (1). Confidence levels can be read from figure^. 



Label 


Po(s) 


P(s/s) 


Ax 2 


(1) 


8.39109273 


5.409 x 10" 14 




(2) 


8.39109148 


3.749 x 10" 14 


1.3 
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Figure 3. The 68% and 90% MLP contours for P and P 
to the complete dataset except for the R98 data, for the 
two solutions (1) and (2). 



We have folded the data on both (P$ , P) solutions (1) 
and (2) in table || to check the relative phasing of all in- 
dividual runs (figure |]). We have then performed the a 



Figure 4- The datasets folded on the (Pq , P) solution (1) 
(left) and (2) (right). 



plication is that RX J0720. 4-3125 is unlikely to be spin- 
ning down under propeller torques ( Alpar 2001 ). In this 
case, under the assumption that the X-ray luminosity of 
the source (Lx ~ 2 x 10 31 di o e rg cm" 2 s" 1 , wher e rfioo = 



d/100pc and d is the distance, Haberl et al. 1997 ) is sup 
plied by energy dissipation, it should be: 

s 



2 x l(n n d?nn < P < 2 x 10 _y rf 



9 j2 



(1) 



The scenario is still consiste nt with the spin-down re - 
cently measured for RBS 1223 ( Hambaryan et al. 2001 ), 
but the value of P reported here for RX J0720.4-3125 is 
well below this range. That also make less plausible an 
interpretation of the hardness ratio profile in terms of ab 



posterivri check with the BuppuStix daLa, which phase on 
correctly — with sululiuii (1), buL iiuL (2) (see the luwesL 

panel of figure [|). This suggest to select solution (1) as 
the most likely timing parameters. In any case, for the 
purposes of our further discussion, the difference between 
the P in (1) and (2) is not significant: both acceptable fits 
to the data have 3 x 10" 14 < P < 6 x 10~ 14 . This is the 
most accurate spin-down measure presented so far for a 
dim NS and, for the first time, it allows a discrimination 
between the proposed models. 

3. Discussion 

The refined value of P reported here is consistent with, 
but two orders of ma gnitude lower than t he extrema of 
the range reported by Haberl et al. (1997) . The first im- 



sorbing matter, co-rotating in the magnetosphere (Crop- 



per et al. 2001 ) . The observed behaviour is more probably 



explained by the angle-dependent properties of the emit- 
ted radiation. 

On the other hand, the slow-down rate of RX J0720.4- 
3125 is still considerable. The other plausible mechanism 
which may account for such large and stable value of P 
is magnetic breaking. For a dipolar magnetic field P ~ 
10" 15 (B/10 12 G) 2 jP s/s, which gives B — 2.13xl0 13 G[] 
A scenario in which this source is powered by accretion 
from the interstellar medium must be therefore ruled out: 
for the present values of P and B the co-rotating magne- 
tosphere will prevent the incoming material to penetrate 
below the Alfven radius. 



1 Here and in the following we specify the discussion to so- 
lution (1) of table ^ 
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B-Decay Mechanism 


Bo 


age 




10 13 G 


(years) 


Hall Cascade 


119.2 


4.5 x 10 4 


Ambipolar diffusion, irrotational mode 


1.9 


3.3 x 10 6 


Ambipolar diffusion, solenoidal mode 


4.1 


1.6 x 10 6 



icantly influenced by decay through ambipolar diffusion 
(Heyl fc Kulkarni 1998|), thus, as in the case of constant 



Table 3. Predicted source age and primordial field for three 
different mechanisms of decay, simulated as in Colpi et al. 
2000. The present values of P and P are those of solution 
(1) in table |^. In all cases, the source is assumed to be 
born with P = 1 ms. 



The corresponding spin-down age is i s( j = P / (2-P) ~ 
2.48 x 10 6 yr, which, given the numerous uncertainties, 
is marginally compatible with that inferred by the cooling 
curves (a few 10 5 yrs for a surface temperature of ~ 80 eV, 
Kapmukci' el al. 20014 iKamiiikei 1 cl al. 20011], |Schaab 



9§?rSchaab ai al. V3§§). The discrepancy in fesg 



e.g _ 
et al 

significant if we notice that what we are probably observ- 
ing in the X-ray is a region of limited size which is kept 
hotter than the average star surface, as inferred by the 



B discussed above, the scenario is compatible with the ob- 
served luminosity. The larger age is also compatible with 
the absence of a remnant. 

If our conclusions are valid, the connection between 
dim INS and AXPs is not so obvious. RX J0720.4-3125 
has a strong, but not superstrong, field which is compat- 
ible with those of the canonical radio-pulsars which have 
passed the death line. On the other hand, having excluded 
accretion, what mechanism causes an X-ray emission con- 
centrated in a fraction of ~ 60% of the star surface remain 
a mystery, as well the related question about the validity 
of using the observed blackbody temperature to locate the 
source in the cooling diagram. The variation of the sur- 
face cooling temperature with the latitude predicted so 
far for stro ng fields ( Greenstcin fe Hartkc 1983] , Possenti 



ct al. 



analysis of the pulse-shape (Cropper et al. 2001 

On the other hand, t s d is representative of the true 
age of the source only in the case in which the magnetic 
field remained almost constant during the star evolution. 
The same condition applies for the validity of the cooling 
curves mentioned above, which do not include the extra 
input of energy released in the neutron star in case of 
field decay. It is therefore of fundamental importance to 
address the field evolution. There are three mechanisms 
which are typically proposed for inducing field-decay: am- 
bipolar diffusion in the solenoidal or irrotational mode 



and H all cascade (Goldrcich fc Reisenegger 1992 , Heyl fc 



Kulkai ju 1998| , JColpi et al. 200C|). By using the simplified 
expressions of Colpi et al. (2000) for the decay laws, we 
have estimated the source age and the value of the mag- 
netic field at the birth of the neutron star, Bo (table ||). 

As we can see, allowing for a mechanism involving very 
fast decay, such as the Hall cascade, we find that the source 
is now ~ 4x 10 4 yr old, and is born with a superstrong field 
Bq ~ 10 15 G. Such a young age is only marginally compat- 
ible with the absence of a remnant and, more important, 
is not compatible with the low observed X-ray luminosity. 



Undcrluminous models have been presented by Kaminkcr 



et al. (2001a), who accounted for the enhanced neutrino 
cooling in presence of strong neutron superfluidity. These 
solutions may match an age of ~ 10 4 yrs for RX J0720.4- 
3125, but, as discussed by the same authors, they must 
probably be rejected since they fail in the comparison with 
observational data of a sample of other neutron stars. 

On the other hand, both mechanisms involving am- 
bipolar diffusion predict a magnetic field quite stable over 
the source lifetime and close to the actual value. Accord- 
ingly, the predicted age is ~ 10 6 years in all cases, close 
to t s d- Below ~ 10 14 G the cooling curves are not signif- 



1996 ) is too smooth to explain the (relatively) small 
size of the emitting region, and different explanations are 
required. 
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